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Abstract 

TITAN (TRIUMF's Ion Traps for Atomic and Nuclear science) is an online 
facility designed to carry out high-precision mass measurements on singly and 
highly charged radioactive ions. The TITAN Penning trap has been built and 

optimized in order to perform such measurements with an accuracy in the 
sub ppb-range. A detailed characterization of the TITAN Penning trap is 
presented and a new compensation method is derived and demonstrated, 
verifying the performance in the range of sub-ppb. 
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1. Introduction 

Penning traps liave proven to be tlie most precise devices for mass spec- 
trometry, both for stable and unstable isotopes [T] . For unstable isotopes with 
half-lives between 5 ms and a few seconds, performing such measurements is 
a challenging task because the unstable isotopes produced via nuclear reac- 
tions need to be delivered in a fast and efficient manner while the subsequent 
measurement need to reach the desired precision. However, the study of sev- 
eral phenomena far from stability benefits from a precise mass determination. 
These physical processes includes: change in the nuclear structure ([2]- [5])) 
the determination of the exact path of the -r, -rp and -z/p processes ([6]- [8]), 
the improvement of the halo nuclei charge radius precision ([9]-[l5]) and the 
testing of the CVC hypothesis ([16] and references therein). The needed rel- 
ative uncertainty on the mass determination for these various cases varies 
from bmjm ~ 10^^ to 10~®. To reach such precision while being accurate, 
great care need to be taken on identifying and minimizing the various sources 
of systematic errors on the mass determination. 

There are currently several experiments dedicated to precise mass mea- 
surements of short-lived nuclei including ISOLTRAP [H] at ISOLDE/CERN, 
CPT HE] at ATLAS/ANL, SHIPTRAP [IH] at GSI, LEBIT [20] at NSCL/MSU, 
JYFLTRAP [21] at JYFL, TRIGA-TRAP at TRIGA Mainz [22] and TITAN 
[23112^ at ISAC/TRIUMF. These experiments are complementary since they 
are set-up at different production facilities and all have a specific reach and 
access to isotopes. TITAN as for example succeeded in measuring several 
masses in the light mass region, including ^He [13], ^Li [2S], ^'^'^^Li [13], 
9'10'iiBe [U] and ^^Be [2^. 

Although the TITAN mass measurement program thrived in measuring 
the masses of very-short-lived (as low as 8.8 ms for ^^Li) halo nuclei, the pos- 
sibility of performing high-precision mass measurements on highly charged 
unstable ions (HCI) is a distinctive feature of TITAN. As it was previously 
demonstrated with stable species at the SMILETRAP experiment [2Z] , HCIs 
are used because the precision of mass measurements performed using Pen- 
ning traps linearly increases with the charge state [2^ . 

Over the past 30 years, extensive work has been done to identify the 
various factors that limits the precision and accuracy of Penning trap mass 
spectrometry. These sources of systematic errors includes: magnetic field 
inhomogeneities, misalignment with the magnetic field, harmonic distortion 
and anharmonicities of the trapping potential, temporal fiuctuations of the 
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magnetic field, relativistic effects and ion-ion interactions. These effects are 
described in the hterature ([28] -[31]) and studied for specific Penning trap 
spectrometers [321 EZl |33] . In this paper we describe how these systematic 
effects where minimized for the TITAN Penning trap and we give an estimate 
for their upper values. Then we determined experimentally the so-called 
mass-dependant systematic shift of the measured frequency ratio which is 
a combination of all the effects that are found to depend linearly on the 
difference in mass-to-charge ratio between the calibrant and ion of interest. 



2. TOF-ICR Penning trap mass measurements basics 

The basic principal behind Penning-trap mass spectrometry [M] consists 
of measuring the cyclotron frequency 

1 qB 

of an ion of mass M and charge g in a magnetic field B. Knowing the field 
strength (this requires a reference mass, as discussed below) and charge state 
of the ion, one can then obtain its mass. In order to reach a high precision, on 
the order of 6m/m < 5xl0~^ on the ion's mass, several requirements need 
to be met. Principally, a magnetic field that is homogenous in the region 
the ion is stored and a sufficiently long observation time are needed. These 
requirements are more easily fulfilled by confining the trapped ions in a small 
volume. 

The Penning trap [35] (schematic shown in figure [T]) is a type of ion trap 
that achieves this using a strong homogenous magnetic field overlaid with a 
quadrupolar electrostatic potential. This potential is created by applying a 
potential difference Vq between a set of electrodes, typically orthogonal hyper- 
boloids of revolution: one forming the ring electrode and the other forming 
the end cap electrodes. The resulting potential in such a configuration is 
given by: 



V2{z,r) = ^y-ry2), (2) 

where {z, r) are the axial and radial coordinates and do = Zq/2 + r^/A is 
typically defined as the characteristic length of the trap. The parameters [zq, 
To) are the distances from the trap centre to the electrodes as defined in fig- 
ure [1} This electrostatic potential configuration axially traps the ions, while 
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Figure 1: (Colour on-line) Schematic diagram of a Penning trap including: the hyper- 
boloids forming the ring electrode (red), the hyperboloid forming the end caps (blue) and 
an axial magnetic field. 

the radial confinement is provided by a magnetic field B parallel to the end 
cap axis. Note that similar confinement properties can also be achieved using 
a cylindrical sets of electrodes that have been orthogonalized [SB] instead of 
the hyperbolical ones. 

There exist analytical solutions for the ion motion in a Penning trap that 
are extensively studied in the literature (see for example [29] and [37|). The 
ion motion is composed of three eigenmotions: one axial of frequency z/^ = 
^y qVo/{Mdl)/{27T) and two radial of frequencies = i'c/2 ± ■y^u'^/A — z/2/2, 
called reduced cyclotron (z/+) and magnetron (z/_). The three eigenfrequen- 
cies typically have the following hierarchy: z/4. ^ z/^ ^ i/_ and it can be 
noted that the cyclotron frequency is not an eigenfrequency, but for an ideal 
(purely quadrupole) trap the relation 



Penning-trap mass spectrometers that uses the time-of-flight ion-cyclotron 
resonance (TOF-ICR) technique |3H] make use of equation ([s]) in order to ob- 
tain the cyclotron frequency of the ion and ultimately, using equation ([T|, 
its mass. It should be stressed that equation ^ is only valid for an ideal 
trap. This paper will present estimates of the various sources of deviations 
from the ideal trap and their effects on the measured cyclotron frequency, 
that ultimately affects the accuracy on a mass measurement in this case, at 



(3) 



holds. 



TITAN. 
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At TITAN, the TOF-ICR technique is used to measure the cyclotron 
frequency. In this technique, the ion's eigenmotions are excited by applying 
a radio-frequency (RF) field. Two types of excitations are typically used: 
dipole and quadrupole. The application of the first one in the radial plane 
at either v_ or z/+ results in the excitation of the corresponding eigenmotion. 
The application of a quadrupolar excitation at the sum frequency z/_ + 
results in a beating between the two different modes. Therefore, by the 
application of an RF excitation at the frequency on an ion initially in a 
pure magnetron motion, a complete conversion of the ion's motion into a pure 
reduced cyclotron motion will occur when i/g = z/_ + A full conversion 
only happens when the RF excitation amplitude Vq and time Tq are related 
by [39]: 

V, = (4) 

where a is the distance from the trap centre at which the RF field amplitude 
equals Vq and ?7o = 1, 3, 5, ... are integer values of the conversion factor r\ 
that allows a full conversion. 

Because ^ z/_, a conversion will result in a gain of the ion's kinetic 
energy E^. in the radial plane. The ion's cyclotron frequency can be derived 
by determining the excitation frequency that yields the largest increase in 
the ion's kinetic energy. This increase is found by releasing the ion from 
the trap and measuring its flight time taken to reach a detector situated 
outside the strong magnetic field region. On the way to the detector, the 
interaction of the ion's motion magnetic dipole moment /i (resulting from 
the RF excitation) with the magnetic field gradient dB^/dz induces a force 
F = V{fl ■ B) = —^^^z that axially accelerates the ions. Since F oc Er, 
the acceleration is the greatest when z/g = i/_ + z/_|_, yielding a shorter time of 
flight. The time of flight as a function of the excitation frequency is described 
analytically by the following integral |39j : 

where zq and zi are the initial and final position of the ion prior and after 
ejection, Eq is its kinetic energy upon leaving the trap, V{z) and B{z) are the 
electrical potential and magnetic fields along the path of the ion. Repeating 
the injection, excitation, extraction, and time-of-flight measurement process 
for different frequencies produces a time-of-flight spectrum such as the one 
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Figure 2: ^Li+ cyclotron frequency resonance taken with a 900 ms excitation time. The 
sohd hue is a fit of the theoretical line shape [39] to the data. 



in figure |2] from wliicli tlie cyclotron frequency is derived as the centroid of 
the minimum. 

To measure the magnetic field of the Penning trap a measurement of the 
cyclotron frequency of calibrant ions has to be carried out. Typically the 
calibrant mass is more precisely known than the mass of the ion of inter- 
est. Hence, the measured quantity from which the mass is computed is the 
frequency ratio 

r, T^c,cal Qcal ' M 

R = = — ^ — • (6) 

Note that the calibrant ion cyclotron frequency Vccaiib. value at the time of the 
z/c measurement is approximated by a linear interpolation of two calibration 
measurements enclosing the measurement of the ion of interest. 

In comparative mass spectrometry, the quantity of interest is the atomic 
mass of a neutral atom, which is given by: 



Qcal 



■ R ■ {meal - Qcai ■ + Be^cal) + Q ■ - 



(7) 



where R is the weighted mean of all measured frequency ratios, B^^cai and 
Be are the calibrant 's and ion of interest's electron binding energies, qcai and 
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q are their respective charge states, mg is the electron mass and rricai is the 
cahbrant atomic mass. The statistical uncertainty on a mass measurement 
is given by the following relation ^01 120] : 

5m 5vc 7 ■ 



m q-B-Tg- 



(8) 



where Nion is the total number of ions detected for the measurement and 7 
is an experiment-dependent constant given by the quality factor [321 120] of 
the time-of-flight resonance spectra. 

Equation ^ describes the precision on a mass measurement using the 
TOF-ICR technique. However, ultimately the accuracy is limited by sys- 
tematic errors that arise from a number of factors. For instance, the trap 
electrodes do not extend to infinity and are truncated. Also, holes in the 
two end-cap electrodes are required to inject and extract the ions from the 
trap and hence disturb the ideal potential. Moreover, the ideal trap assumes 
perfect geometrical alignment of all applied electrostatic and magnetic fields. 
In reality, misalignments between each trap electrode and distortion in the 
shape of the electrodes exist due to technical limitations in the achievable 
machining tolerances and affect the trapping potential [29l [30]. There are 
also misalignments of the trap's electrode structure principal axis with the 
magnetic field axis, and deviations of the magnetic field in the trapping re- 
gion. Other effects arise due to the Coulomb interaction between stored 
ions. Moreover, fluctuations of the magnetic field strength over time, and 
relativistic effects have to be taken into account. These various effects result 
in a different measured cyclotron frequency from the true Vc given by equa- 
tion ([1]). The resulting frequency shift modifies the measured frequency ratio 

■Rmeas. • 

jj ^c,cal ~l~ ^^c,cal /q\ 

from the ideal frequency ratio Rideai = ^c^cail'^c- 

The large value of the cyclotron frequency, in the MHz range, compared to 
the frequency shifts Az/c, in the Hz range, allows one to state that Az/c/z/c <^ 
1. This consequently leads to a relative frequency ratio shift of 

A-R Rmeas. Rideal ^^c,cal Az^c (10) 

R Rideal ^c,cal. 

For most systematic effects studied in this paper, Auc^i ~ Az/c 2 = Az/^, there- 
fore the relative frequency ratio shift will typically have the form AR/R = 
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Figure 3: (colour on-line) The TITAN experimental setup which includes an RFQ, a high- 
precision Penning trap, an EBIT, a time-of-flight gate and an off-line ion source, a) Shown 
in solid-red is the path of the beam when mass measurement on singly charged ions (SCI) 
is performed, b) In dashed-blue is the path for mass measurements on highly charged ions 
(HCI). 



(27r ■ Az/c/-B) ■ A(m/g), where we defined A(m/g) := rricai./ qcai. — From 
equation ( [Io| ), two main conclusions are drawn. Firstly, relative frequency 
ratio shift are in general smaller than the individual relative frequency shifts. 
Secondly, by measuring the frequency ratio of two species of similar mass-to- 
charge ratio m/g, one can reduce AR/R. 

3. The TITAN mass measurement Penning trap 

The high-precision mass measurements carried out at TITAN (shown 
in figure [s]) are achieved through a series of steps. First, the continuous 
ion beam from ISAC (Isotope Separator and Accelerator) |1T] is delivered 
to TITAN where it is cooled and bunched using a gas-filled linear radio- 
frequency quadrupolar (RFQ) trap [12] • The subsequent step varies depend- 
ing on whether a mass measurement is performed using singly charged ions 
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Figure 4: (colour on-line) Schematic of the TITAN Penning trap electrode configuration 
formed by the hyperbolic ring (labeled (1)), end cap electrodes (2), tube (3) and guard 
(4) correction electrodes. The RF is apphed on (4) and the blue-red color code express 
the opposite phases of a quadrupolar excitation. The characteristic dimensions are given 
m table [H 



(SCI), or highly charged ions. The choice depends on the required preci- 
sion, half-hfe and production yield. The ions can either be transferred to an 
electron-beam ion trap (EBIT) (|13]-|15]), where charge breeding takes place 
(blue path in figure [3]), or sent directly to the Penning trap (MPET) where 
the mass of the ion of interest is determined (red path in figure [3]) . 

The TITAN Penning trap is shown in figure |4| with characteristic dimen- 
sions given in table [Tj The trap is composed of two hyperboloids of revolution 
forming one ring (label (1) in figure |4]) and two so-called end-cap electrodes 
(2). The ions are axially trapped by a harmonic quadrupole electrostatic po- 
tential produced by a potential difference, Vq, between the ring and the end 
cap electrodes, as shown in figure |4j Some anharmonicities in the trapping 
potential are introduced by the holes in the end-cap electrodes and by the 
finite size of the hyperbolic electrodes. Two sets of correction electrodes (la- 
beled (3) and (4) in figure |4]), are used to compensate for higher-order electric 
field components (for more detail see section 3.3). The radial confinement is 
provided by a magnetic field B. This section describes how various sources 
of systematic errors effects are studied and minimized. 
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Table 1: Characteristic dimensions of the TITAN Penning trap shown in figure [4j 



Dimension 


Value (mm) 


'"o 


15 


zo 


11.785 


do 


11.21 


H 


41.48 


D 


61 


d 


4 



3.1. Spatial magnetic field inhomogeneities 

For an ideal Penning trap, one assumes that the magnetic field strength 
is constant across the trapping region, i.e. B{x,y,z) = Bq. For real traps, 
however, magnetic field inhomogeneities are created by the finite size of the 
solenoid and magnetic field distortion due to the magnetic susceptibilities of 
the trap material [30] or imperfections in the solenoid due to the finite size of 
the coil wire. In all cases, the lowest-order contribution to the magnetic field 
axial projection inhomogeneities has a quadrupole component, as seen in 

B,{z,r) = Bo{l + P2{z'-ry2)} (11) 

where z is the ion axial oscillation amplitude, r is the ion radial position, 
Bq is the unperturbed magnetic field strength and is the strength of the 
lowest order inhomogeneity component. This constant has been emperically 
determined for various other systems and typically ranges from 10^^° to 10~^ 
mm-2 for TOF-ICR Penning traps |30] and [20]. 

The specifications for the construction of the TITAN magnet were such 
that the Bq = 3.7 T field is homogenous within AB^/Bq = 1 part-per-million 
(ppm) inside a 2 cm long by 1 cm diameter cylinder about the trap centre. 
Also, the trap electrode structure was manufactured by minimizing the re- 
quired amount of material to realize the geometry. Moreover, only material 
of low magnetic susceptibility such as high conductivity oxygen-free copper 
for the electrodes and sapphire for the insulators was used. 

With these considerations, one can estimate the effect of the magnetic 
field inhomogeneities on the measured frequency ratio. The frequency shift 



from the magnetic field inhomogeneity (11) is [30 



Az., = (32U^ I {z' - r^) - ^ {rl + r!) | , (12) 
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where r_|_ and r_ are respectively the ion reduced cyclotron and magnetron 
radii. Because z/c ^ z/_, the term in braces is weakly mass dependent, leading 
to an overall effect on the cyclotron frequency ratio of: 

(^) --/3.K + r^)-^A(m/g), (13) 

where we used i'_ ^ Vo/(47ri?(io). The size of the TITAN Penning trap end 
cap hole can be used as estimate of the maximum workable or r^, from 
which we estimate (r^ + r^) < 4 mm^. Using this estimate and assuming 
= 1 X 10~® mm~^, one gets an absolute value in the shift in the frequency 
ratio of (AR/R) mag.inhom. < 1-2 X 10"" -Vo ■ A(m/g), where A(m/g) 
is in units of amu/charge. Assuming Vq = 35.7 V, the uncertainty due to 
the spatial magnetic field inhomogeneities becomes {AR/ R) mag.inhom. < 4.3 x 
10^^°- A(m/g), which is one order of magnitude below the standard precision 
aimed for the TITAN Penning trap system. 

3.2. Harmonic distortion and misalignment of the magnetic field axis 

The ideal Penning trap assumes a perfect alignment between the trap 
electrode structure axis and magnetic field axis (i.e. B = Bqz). It also 
assumes that the electrodes are aligned with respect to each another and 
without surface imperfections. In reality (figure |5] (a)), the magnetic field 
could have some misalignment with the trap axis. Also, the trap electrodes 
could have deformations as shown in figure |5] (b), and be misaligned with 
respect to one another. This would lead to a finite asymmetry parameter e, 
resulting in a distorted potential [28] given by 

Vharm.dist. = ^ {(^ + e)^' + ' ^)y^} ' (l^) 

These two imperfections modify the equation of motion of the ion in the 
trap resulting in a change of their eigen frequencies, modifying the measured 
cyclotron frequency according to [28] : 

A..= (V-ie^).z7_, (15) 

where Auc is the cyclotron frequency shift, T'_ is the measured magnetron 
frequency and 6 is the angle between the trap and magnetic field axis. The 
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Figure 5: a) Schematic of the electrode structure ahgnment with the magnetic field axis, 
b) Top view of a elliptically distorted ring electrode that leads to a non-zero asymmetry 
parameter e. Also shown is the undistorted ring (dashed lines). 



corresponding frequency ratio shift is given by: 

(Afi/fi)„„,^(V-l.^).(^).(Jl-), (16) 

where z7+,cai is the measured reduced cyclotron frequency of the cahbrant. 



By approximating v+^cai ~ T-'c,caU equation (16) can be simphfied to 

m/RUs. = i^-^e' - \e'^ ■ ^ ■ A(m/g). (17) 

The angle of misalignment 6 was minimized by a precise alignment of the 
Penning trap electrode structure with the magnetic field axis using an elec- 
tron beam. In addition, the system was built by requiring tight machining 
tolerances of lO/xm for the trap electrodes and insulators. The trap vacuum 
chamber has been carefully aligned using an electron source positioned at 
the trap centre location, a phosphor screen placed at the end of the vacuum 
chamber, and three concentric targets as shown in figure |6j The targets were 
made of an aperture with two metal stripes accurately positioned to form a 
cross within 0.01 mm of the hole centre. The stripes were 0.1 mm thick and 
the hole was 8 mm in diameter. The three targets were then secured in place 
along the Penning trap optics support frame, which in turn is centred along 
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shadowed "x"s strips forming an "x" 

inside a co-centric hole 

Figure 6: Schematic of the Penning trap vacuum chamber ahgmncnt using an electron 
source. 



the vacuum tube ensuring the co-centricity of the targets. The vacuum tube 
is made using a pulled honed titanium tube with a very low tolerance on 
centricity of less than 0.13 mm over the length of 1.23 m. The tube acts as 
an optical bench for the trap structure. The cross patterns were rotated from 
one target to another as shown in figure |6] in order to facilitate the alignment 
of the chamber. 

The alignment principle is based upon the fact that in the magnetic field, 
the electrons are guided along the field lines. A proper alignment of the 
vacuum chamber co-centric with the magnetic field lines is ensured when the 
three shadowed images of the targets metal strips are aligned and when a 
circular image of the electron beam spot is observed on the phosphor screen. 
Therefore, once the chamber has been adjusted such that the electron beam 
can pass through the three apertures, which are placed along the tube axis, a 
fine alignment is performed by aligning the "shadows" of the three crosses on 
a phosphor screen. This is done by moving the chamber with respect to the 
magnet housing in the x-y direction using a fine-thread external mechanical 
alignment mechanism. 

The chamber has been finely adjusted until a displacement not worse 
than 0.04 mm between the three targets has been reached. Considering that 
the distance between aperture one and three (see figure [6]) is 590.5(1) mm, 
this gives an upper limit on the misalignment of the vacuum chamber with 
respect to the magnetic field axis of Ochamber < 7 x 10^^ rad. 

An additional source of error comes from a misalignment of the trap 
electrode structure with respect to the support frame. The Penning trap 
electrodes are isolated from one another and from the support frame using 
sapphire spheres that are placed in a depression of the electrode. The entire 
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Figure 7: a) Left: one-piece ring electrode. Right: gold-plated Penning trap electrodes 
including the sapphire spheres (shown on top of the sliced guard electrode). Shown is a 
Canadian dollar coin for scaling, b) The TITAN mass measurement Penning trap (MPET) 
placed in the support structure frame. The trap structure is held in place by compression 
using two support rings. 




Figure 8: Schematic to demonstrate misalignment stemming from unequal sapphire sphere 
and holding hole size leading to a misalignment between the trap and magnetic field axes. 
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structure is held in place by compression (see figure [7j), and the support 
structure is connected to the vacuum tube. The maximal misalignment of 
the trap with respect to the support structure happens when the sapphire 
balls are too wide, while the holes in which they are located are too small. 
This combination with extreme opposites for the mechanical tolerances on 
either side of the electrode structure, as shown in figure |8| would lead to a 
maximal shift. 

We consider sphericity tolerances on the sapphire balls of 5 /im and toler- 
ances of 20 /zm for the grove in which the balls are placed, which correspond 
to a doubling of the specified manufacturing tolerances. In the worst case, 
both the correction tube electrode with respect to the support frame and 
the end cap electrode with respect to the tube electrode are misaligned. Un- 
der these conditions, the misalignment is 5 = 90 fim. Considering that the 
support sphere are placed along a circle of radius i? = 21.6 mm, the largest 
possible tilt is 

Osupp.ma.. = 5/R = 4.2 X 10-3. (18) 

Since the error stemming from the machining tolerances dominates [Osupp. S> 
(^chamber) the error from the chamber alignment, the upper value on 6^, is 9max 
= 4.2 X 10-3. 

A non-zero asymmetry parameter, e, is caused by a number of effects. 
We discuss them and show how they can be minimized. Firstly, localized 
oxidation patches on the surface of the electrodes cause undesired inhomoge- 
nous electric fields because of the different dielectric constant of the patch 
[16] . Such effects are minimized by gold-plating the trap electrode surfaces 
(see figure [T]). Secondly, the quadrupole deformation of the electric poten- 
tial in the xy-plane is minimized by applying the RF-field on the correction 
guard electrodes to avoid splitting the ring electrode. Thirdly, misalignment 
of the ring electrode with respect to the trap axis is minimized by using high- 
tolerance sapphire spheres on which the trap electrodes sit and by requiring 
and ensuring tight electrode machining tolerance. 

Based on the manufacturing specifications, and the required machining 
tolerances, the maximum value for the ring electrode tilting angle a is found 
to be 0.0016 radians. Applying trigonometry and the equation for the ring 
electrode hyperbola, the asymmetry parameter e due to tilted ring electrode 
is 

euit = sin^ a ■ (1 + (ro/^o)') = 1-1 x 10"^ (19) 
Lastly, elliptical deformation of the Penning trap electrodes would cause 
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a non-zero asymmetry parameter e. Elliptical deformation of the ring elec- 
trode will have the largest impact on the electrostatic potential, as it is the 
closest electrode to the trap centre. An elliptical deformation of the ring 
electrode (see figure |5|^b)) corresponds to an elliptical Penning trap, a spe- 
cial case of hyperbolical Penning trap that has been studied extensively both 
theoretically |17j and experimentally [IH]. The elliptical deformation of the 
ring electrode can be described by: 

+ . = 1, (20) 



e 



where e gives an ellipticity that varies from < e < 1. Assuming the ring 
electrode radius from the trap centre is tq = 15 mm, and conservatively 
doubling the machining tolerance 5 = 0.01 mm, we get 

emax = — = 2.6 X 10-^ 21 

The error on the frequency ratio due to both the asymmetry parameter e and 



the angle 9 is given by equation (17). The maximal error on the frequency 
ratio is obtained when e = and 6* = 4 x 10"'^ and is equal to {AR/ R)mis. = 
1.2 X 10^^°V~^ ■ Vq ■ A(m/g). Using Vq = 35.7V, one obtain an upper value 
of {AR/R)mis. = 4.3 X lO-^V"^ • A{m/q). 

3.3. Non-harmonic imperfections of the trapping potential 

The holes in the end cap electrodes and the truncation of the Penning trap 
hyperboloid structure cause the trapping potential to be non-ideal and hence 
non-harmonic. Therefore one needs to consider octupole and dodecapole 
corrections to the trapping potential [30]. These corrections are given by: 

^0 \ I 6 15 4 2 , 45 2 d 5 g 



where C4 and Cq are the octupole and dodecapole correction strengths, re- 
spectively. Because of the dependence of the potential, the contribution 
of higher-order terms become increasingly smaller. The procedure to calcu- 
late the frequency shifts due to the non-harmonicities is given in detail in 
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[22] and [30], and result in shifts in the radial eigenfrequencies of: 

Az/± ^ ±7§^- {{rl + 2rl) - 2z'] (24) 
Ai/± ^ i^^''- + ^'^'(^i + - (4 + 3^ + 6r^r2 )} , (25) 

for the octupole and dodecapole terms, where the frequency shifts are nearly 
mass independent. This is because the magnetron frequency is very weakly 
mass-dependent [39] . Substituting z/c = + u_ gives: 

A.. ^ \^^^^^--{C. + \^ (3.^ - - r^)}. (26) 

Assuming that no compensation voltage is applied on the correction-tube 
and -guard electrodes (see figure |8]), one obtains the following coefficients: 
C4 = 0.004 and Cg = -0.082 for the TITAN Penning trap. Using a repre- 
sentative oscillation amplitude of 2; = 3 mm and (r^ + r^) = 4 mm^, we 
get AR/Rpot.inhom. = 2.8 X 10~^ • Vq ■ A(m/g). This is over 20 times larger 
than any other previously discussed frequency shift and for a typical trapping 
potential of Vq = 35.7 V, this can lead to a frequency ratio shift of lxlO~^. 
Therefore, in order to perform accurate mass measurements at the level of 
Sm/m ^ 5x 10~^, it is necessary to minimize the non-harmonic coefficients. 

3.4- Additional sources of systematic errors 

Other sources of systematic errors include magnetic field fluctuations in 
time, ion-ion interaction and relativistic effects. These effects have been dis- 
cussed for the speciflc case of the ^Li mass measurement (see f25[ HH] for 
a detailed description of these effects). The systematic error on the lin- 
ear interpolation of the calibrant cyclotron frequency due to time-dependant 
fluctuations of the magnetic fleld was found to be 0.04(11) ppb/h [25]. In 
most cases, the calibrations are spaced by less than one hour. Therefore, this 
represent a small systematic error. 

Another source of systematic error would come from fluctuations of the 
trapping potential over time. This effect was investigated by monitoring at 
the potential applied to the speciflc electrodes at two occasion over a two 
months interval. The largest observed change in potential were for the ring 
electrode: Akring = 0.0022(1) and the injection-side correction tube electrode 
Aktube = -0.0009(1). Such change in ktube from the optimal ktube = 1.53(2), 
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would results in frequency shift between rj = 0.5 and 1.5 of 3.4(4) mHz 
for ^Li^, which is much lower than the 80(50) mHz change observed in the 
previous section. Therefore, trapping potential fluctuations are considered 
to have a small relative effect. 

4. Compensation of the Penning trap electrical potential 

From the previous sections, the largest possible source of error on the 
measured frequency ratio would come from the non-harmonic terms in the 
trapping potential. Therefore, these terms need to be minimized as they can 
induce a large shift in the cyclotron frequency. This is achieved using the 
correction guard and tube electrodes shown in figure |4} 

The optimal correction guard and tube potentials with minimal non- 
harmonic coefficients Cn>2 have been estimated through the chi-square min- 
imization of the difference between the potential produced by the trap elec- 
trodes and a quadratic target potential: 



X 



2l2 



$^{K.(.)-(^Ao)^} 



(27) 



The effective axial potential Vax is a linear combination of the axial potentials 
produced by the individual electrodes: 



^cap^capi^z) ~l~ ^ring^ringi^z) ~l~ ^guardVguardi^^) ~l~ ^tubeS^tubei^^^ i (^8) 



where Vi{z) is the axial potential produced when IV is applied on the surface 
of a given electrode and fcj are the scaling coefficient determined by the 
chi-square minimization. The potentials Vi{z) corresponding to the TITAN 
Penning trap geometry were obtained using the Laplace equation solving 
capabilities of the ion-optics simulation software SIMION [50] . Upon solving 



Table 2: Calculated normalized potential ki needed to be applied on the ring, correction 
tube and guard electrodes in order to optimally compensate over a range of 8 mm from 
the trap centre. 



h ■ 

luring 


ktube 


^guard 






-0.786 


1.640 


0.078 


-7x10-6 


5x10-5 



equation (27) with k, 



^cap 



1, one obtains the optimal potentials shown in 



18 



table |2} The size of the residual C4 and Cq coefficients for the optimal 
potential configuration is obtained by the least-square regression of 

^0 / ^ , ^^2 2 I ^4 4 Ce 6 



The resulting coefficients shown in table |2] are a factor 1000 smaller than the 
values presented earlier, and would result in a significant reduction of the 
cyclotron frequency shift. 

4-1. Penning trap compensation using a dipole excitation 

The usual procedure [SI] to compensate the electrostatic potential of a 
Penning trap consists of measuring the reduced cyclotron frequency of the 
ion in the trap for two values of the axial oscillation amplitudes z: one with 
z ~ and the other with z > 0. The optimal compensation is the one that 
minimizes the difference between these two reduced cyclotron frequencies: 

5u^ = Uj^{z = zq) — z/_|_(2 ~ 0). (30) 

Since the ions are trapped dynamically, the amplitude of the oscillations 
can be controlled by the closing time of the trap. Assuming the correct en- 
ergy, the ions will have their minimal kinetic energy once they reach the trap 
centre and if at this point the trap is closed, the axial oscillation amplitude 
of the ions should be minimized. However, if the trap is closed at different 
times, earlier or later, it results in larger axial oscillation amplitudes. 



From equations \2A\ and (25), the reduced cyclotron frequency changes 
in a quadratic form with an extremum at the trap centre. This is shown 
schematically in figure [9] (left) for different correction tube potentials. Note 
that the extremum in frequency found for capture time tcap = 40.3 /iS (the 
capture timing is calculated from the extraction of the bunches from the 
RFQ) correspond to the time for which minimal ion axial oscillation ampli- 
tudes are observed. By changing either the correction guard or tube poten- 



tial, one changes the values of the coefficients as expressed in equation (33 ). 



This changes the amplitude and direction of the concavity of equations (24) 



and (|26]) (figure [9] (left)). 

Thus, the Cj>2 coefficients are minimized by changing the potential on 
the correction electrodes and taking the difference between the measured 
Z/+ at the trap centre and the at a location away from the centre, i.e. 

5u+ = V+itcap = 39.8/is) - U+{tcap = 40.3/is). 
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Figure 9: Left: Variation in the reduced cyclotron frequency as function of the capture 
time (hence position of the ions at the potential switching point) in the Penning trap 
for different correction tube settings ktube for ^Li+. Note that close to ktube — 1-40, the 
variation in j/+ seems to be minimal. For these scans, kguard — -0.05. Right: Change in 
the reduced cyclotron frequency difference of ®Li+ with the correction tube potential 
ktube for a correction guard potential of kguard — -0.05. Note the linear change in 5vj^ 
with tcap- The linear fit crosses 5v+ = for ktube = 1.464(4). 



Figure [9] (right) shows the hnear behaviour in the reduced cyclotron fre- 
quency difference 5^+ with the correction tube voltage for kguard = -0.05. A 
linear regression of the data shows that Su^ crosses zero for ktube = 1.464(4). 
This corresponds to one of the possible compensations. Therefore, by repeat- 
ing the procedure for different correction guard potentials kguard one finds a 
family of different compensation settings. 

4-2. Motivation for a compensation using two methods 

In this section we demonstrate that using only one electrostatic potential 
compensation method (such as the one presented in the previous section) 
leads to ambiguous values for the optimal correction tube and guard volt- 
ages, motivating the need for compensating the trapping potential using two 
different methods. 

In order to find the behaviour of the reduced cyclotron frequency differ- 
ence (5z/+ with the correction tube and guard potential, we investigated the 
C4 and Ce coefficients behaviour as function of ktube and kguard- This was 
done by varying the scaling coefficients over the ranges -1.0 < kguard < 1-0 
and 0.8 < ktube < 2.0 and by calculating the C4 and Cq coefficients using 



equation (29). Figure 10 shows that the C4 and Cg coefficients vary linearly 
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Figure 10: Linear variation in the strength of the C4 and Cg coefficients as function of 
both ktube and kguard- Note that the planes cross at kmbe — 1-64 and kguard — 0.08. 



with both ktube and kguard, giving rise to planar surfaces of equations 

C4 = 0.004 - O.OOOSktube - Omikguard (31) 

Ce = -0.083 + 0.050ktube + O-OUkguard (32) 

in the ktube-kguard space. As discussed in [29], the octupolar term C4 is 
mainly corrected by placing a correction guard electrode between the ring 
and end cap electrodes. This is confirmed by the strong dependance of C4 
with kguard showu in (31). Figure 10 shows that the dodecapole term Cg is 



mainly affected by ktube, confirming the literature PU]. 



By inspecting equation (28) and equation (29), it can be shown that the 



linear behaviour of C4 and Cq can also be generalized for higher order Ci. 

Ci = CLiktube + b^kguard "I" Q (33) 

which allows one to write the reduced cyclotron frequency difference as: 



Su+ = ^ a2ih2iktube + ^ h2ih2ik 

guard 



(34) 



1=2 
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where are functions of the axial and radial positions of the ions in the trap 
that are not affected by variations of C4 and Cq. The optimal compensation 
condition = leads to optimal values for ktube and kguard, lying along a 
line given by the equation 

EOO 7 7 ^^00 7 

l^tube sr^oo 7 '^guard -^-^oo 7 ■ v^"^) 

Since there can only be one sets of kmbe and kg^ard that leads to a mini- 



mal value of the C4 and Cq coefficients (see figure 10), and to an optimal 
compensation, one needs a second compensation approach that selects the 
correct setting along this line. This is achieved by carrying out two indepen- 



dant methods of compensating the trap. It should be noted that figure 10 
does not include higher order terms which are present when C4 and Cg are 
effectively zero. However, the terms with the dominant contribution are 
minimized. 

4-3. Penning trap compensation using a quadrupole excitation 

The compensation using a quadrupole excitation consists of measuring 
the cyclotron frequency of the ion for two different conversion factors, rj. 
This factor defines the sizes of the magnetron and reduced cyclotron radii 



(r_ and r+) at the end of the excitation phase and from equation (26), 
changing these radii leads to different cyclotron frequencies. In this method 
the non-harmonic terms are minimized by finding the potential kg^ard and 
ktube that minimizes the change in the cyclotron frequency with rj. The 
conversion factor itself is changed through a variation of the RF amplitude 
Vq for constant excitation time Tg (see equation Q). 

The change of i>c with Vg was studied by numerically solving the equation 
of motion with an added C4 term. Note that the cyclotron frequency is also 
be modified by Cq and higher order terms, but to simplify we only study the 
changes involving C4. Upon solving these equations of motion, one obtains 
the radial energy profile as a function of the detuning frequency Aug = Vg — Vc- 

EriAug) = {xiAugf + y(Az/,)2) , (36) 
where the dots denotes temporal derivatives and (x, y) is the position of 



the ion in the radial plane. Figure 11 shows that for C4 7^ 0, the radial 



energy profile is no longer symmetric. The deformation is more pronounced 
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^^fif ^RF 



Figure 11: Radial energy gain (calculated from equation (36)) as a function of the 
detuning Ai/npTup for conversion factor rj — 0.5 and 1.5 for CiTup = 0.01. Note that 
this is a large C4 value used to demonstrate the effect. The energy gain was normalized 
to the radial energy obtained from C4 — 0. 



for over-converted resonances (17 > 1), due to their smaller line width. For 
the under-converted case {rj < 1), the centre frequency is more shifted. This 
is due to the ion magnetron motion being not fully converted into reduced 
cyclotron motion. Hence the ion then spent more time in regions where the 
Ci>2 components are larger, leading to a larger shift in the centroid frequency. 

Figure 12 (left) shows how the cyclotron frequency Vc changes with rj for 
three different non-zero C4 = -0.003, 0.001 and 0.005. As expected from 
equation (26), the larger C4, the more sensitive with rj the cyclotron fre- 
quency becomes. Also, the shift in frequency Az/c flips sign together with 
C4. Figure 12 (right) gives a similar view of this phenomenon, except C4 is 
varied for r/ = 0.5 and rj = 1.5. This figure shows the only case where the 
cyclotron frequency is the same for the two conversion factor when C4 = 0. It 
also shows that for values of C4 close to zero, the cyclotron frequency change 
linearly with C4 for both rj = 0.5 and f] = 1.5. Therefore, the optimal trap 
compensation that minimizes C4 and other higher order terms will be equal 
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Figure 12: Left: Shift in the cyclotron frequency with the conversion factor rj for C^Trp 
= -0.003, 0.001 and 0.005. Right: Change in the frequency shift Ai^c as function of the 
octupole coefficient C4 with 77 = 0.5, 1.5. These figures are numerical calculations results. 



to the correction guard kguard and tube kmbe setting for which the cyclotron 
frequency difference for rj = 0.5 and 1] = 1.5 equals zero. 

In the following, this compensation is performed by scanning the correc- 
tion tube potential for three different correction guard voltages with values: 
kguard = -0.05, 0.06 and 0.08. For each ktube, the cyclotron frequency was 
determined for rj = 0.5 and r] = 1.5 by fitting the centroid using a Gaus- 
sian curve as shown in figure 13 The previous numerical calculations results 



shown in figure 11, reveal that in the presence of a large C4 term in the 



electrostatic potential, the radial energy gain profile of the ion when i] ~ 1.5 
becomes asymmetric about the cyclotron frequency. Figure 13 (right) shows 
that such effect gets translated into an asymmetric TOF resonance spectrum. 
When the compensation gets better, the TOF resonance spectrum when rj ~ 
1.5 becomes more symmetric as shown in figure [13] (left). A non-optimal 
compensation also results in a large change in the cyclotron frequency for 



7] ~ 0.5, as shown in figure 13 (right). 

Next we calculated the cyclotron frequency difference 6i>c = I'dv = 0.5) — 
i>c{ri = 1.5), from which = was found from linear regression. The 
difference between the under-converted cyclotron frequency and the optimal 
over-converted cyclotron frequency for different correction tube voltage and 
using kguard = 0.06 is shown in figure 14 The linear regression of 6i>c for 



this kguard yielded an optimal correction tube voltage of ktube = 1.53(1). 
Note that the two time-of-flight profiles at ktube = 1-52 and 1.54 for r] = 1.5 
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Figure 13: Over- (77 ^ 1.5) and under- (?; ^ 0.5) converted time of flight resonant spectra 
for kguard = "0.05. Left: A close to optimal setting {ktube — 1-55) resulting in a symmetric 
over- converted TOF spectra. Right: A clearly non-optimal setting {ktube = 1-75) resulting 
in a crocked over-converted TOF spectra. 



were the most symmetric, which is a clear indication of reaching the optimal 
compensation. 

The same procedure was repeated for correction guard voltages of kguard 
= 0.08 and -0.05. The three different optimal kguard and ktube found using this 



method, together with a linear regression, are presented in figure 14 (right) 



This graph also shows the results from the compensation using the dipole 
resonances (see section 4.1). The optimal compensation of the non-harmonic 
terms in the trapping potential was taken as the intersect of the two lines. 
These two lines meet for ktube = 1-53(2) and kguard = 0.01(2). For the typical 
end cap potential Vcap = 20V used at TITAN, these corresponds to Vtube = 
30.6(4)V and Vguard = 0.2(4)V. 

The validity of this result was verified by investigating potential changes 
in the cyclotron frequency with the conversion factor r]. Figure 15 shows that 
the changes in the cyclotron frequency is minimal for ktube = 1-54 and kguard = 
0.02 compensation where the largest change in frequency with the conversion 
factor is 80(50) mHz. For comparison purposes, three other settings are 
shown. When ktube is very different from the optimal value, it induces a 
strong shift in the cyclotron frequency as the conversion factor is varied from 
about 0.3 to 1.6. As expected from figure 12, this shift change direction for 
ktube values above or below the optimal one. 

Figure 15 also shows the setting ktube = 1-46 and kguard = -0.05 which was 
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Figure 14: Left: Change in the fitted cyclotron frequency difference I'dv ~ 0-5) — i^dv = 
1.5) with the scaled correction tube potential ktube for kg^ard — 0.08. Right: Optimal 
kguard and ktube found using a minimization of Sv^ and Svc- The intersection of the solid 
and dashed lines correspond to the optimal combination of kguard and ktube- 

previously found as being optimal using the minimization of (5z/+ method. 
However, when varying the conversion factor t], changes in the cyclotron 
frequency of 260(80) mHz are observed, which are three times larger than 
the ktube = 1-54 and kguard = 0.02 compensation. This observation confirms 
the relevance of compensating the non-harmonic terms of the potential by 
using two different observables. 

5. Experimental determination of the mass-dependent frequency 
ratio shift 

The various systematic effects studied and minimized in section 3.1 to 
3.3 all results in relative changes in the frequency ratio of the form AR/R = 
{2ti ■ Auc/ B) ■ A{m/q). At the end of these sections, we presented upper limit 
estimates of AR/R due to the magnetic field inhomogeneities, the misalign- 
ment of the trap electrodes with the magnetic field, the harmonic distortion 
and the non-harmonic terms in the trapping potential. These estimates were 
based on the chosen trap geometry, machining tolerances and the trap align- 
ment with the magnetic field. 

A more realistic value for the total contribution of the systematic effects 
for which the relative changes in the frequency ratio depends on A{m/q) 
can be evaluated experimentally by measuring the frequency ratio of ions 
with different mass-to-charge ratio. Using these frequency ratio, the atomic 
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Figure 15: Fitted cyclotron frequency as a function of the conversion factor t] for different 
correction tube (ktube) and guard (kg^ard) potentials. 



mass of one of the two species is calculated using equation ([T]) and compared 
to the most precise value from literature. The difference in mass Am = 
m(TITAN) — m(literature) is then used to compute a combined systematic 
shift: 

— ■ ^ = ^ . 1 (37) 

R A{m/q) m A{m/q) 

This so-called mass-dependent shift has been calculated from four differ- 
ent mass measurements: ^Li [l9], ^^Na, ^^K, and ^^K using as calibrant ^Li, 
H3O, ^^Na, and H3O, respectively. These measurements were taken using 
the TITAN off-line ion source and a trapping potential of Vq = 35.7V. The 
resulting frequency ratios are presented in table [3} while the mass difference 
for ^Li, ^^Na, ^^K and ^^K compared to the 2003 Atomic Mass Evaluation 
(AME03) [52j values and a more precise measurement from the Florida State 



University (FSU) Penning trap [HSj are shown in figure 16 All masses agree 
with the FSU measurement within one a. The total shift in the frequency 
ratio was taken as the weighted mean of the three measurements yielding 
-0.2(2) ppb/u. By dividing this value by the trapping voltage used: Vq = 
35.7V, one obtains a relative change in the cyclotron frequency ratio due to 
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Table 3: Frequency ratios R — v^^cailvc of ®Li [55], ^'^Na, "^^K, and using as calibrant 
^Li, H3O, ^^Na, and H3O, respectively. Difference in the mass excess Am of ^Li, ^'^Na, 
and ■^^K as measured with the TITAN and from the FSU [SS* Penning traps. From 
these measurements the systematic error due to the incomplete compensation i^v jiyA) is 
derived. N is the number of measurements taken. 



Specie 




/\{m/q) 


Am (eV) 


syst. error (ppb/u) 


N 




857 332.053 6(37) 


1 


1(24) 


-0.1(2.2) 


43 




1 208 823.886(4) 


4 


-66(77) 


-0.8(9) 


12 


39K 


1 694 844 656(12) 


16 


-52(263) 


-0.1(5) 


13 




2 153 834.195(13) 


22 


-109(226) 


-0.1(3) 


15 


Total 








-0.2(2) 


40 



the total combined systematic effects equal to 

{^R/R)totai = -4(6) X 10-12 . ■ A{m/q). (38) 

Note that if the less precise AME03 masses are used instead, the mass- 
dependant shift becomes -0.6(3) ppb/u and {/S.R/ R) total = —1-7(8) x 10-^^ • 
V^o-A(m/g)). 

6. Summary and outlook 

Penning trap mass measurements to a level of precision and accuracy of 
5m /m ~10-^ are only made possible if detailed systematic studies of the sys- 
tem is performed. Such studies of the TITAN Penning trap are presented, 
in particular the different sources of systematic errors on the measured cy- 
clotron frequency arising from the imperfections of the Penning trap, such as 
the magnetic field inhomogeneities, the misalignment of the trap electrodes 
with the magnetic field, the harmonic distortion of the trap potential, and 
the non-harmonic terms in the trapping potential. 

The total systematic error on frequency ratio determination at TITAN 
was found to be {AR/R)totai = —4(6) x lO'^^ ■ A{m/q) ■ Vo- These mass- 
dependant systematic error depend on the mass-to-charge ratio difference 
A(m/g). This means that when the mass measurement is performed us- 
ing a calibrant and species of similar mass-to-charge ratio, the shift on the 
frequency ratio will be effectively quenched. Also, these estimates are pro- 
portional to the trapping potential Vq, and their contribution to the mass 
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Figure 16: Masses of ^Li [49], ^^Na, ^^K, and measured by the TITAN Penning trap 
using ^Li, H3O, ^'^Na and H3O respectively as calibrant, compared to both the AME03 
[52] and a more recent measurement from the FSU Penning trap |53| . 



measurement systematic error can be reduced by using a small trapping po- 
tential, hence operating in a so-called shallow trap. 

We also presented compensation of the trapping potential using a new 
general method of compensation. We performed this compensation using 
two different methods in order to optimize the combination of correction tube 
and guard voltage that provide the maximal compensation of the trapping 
potential. Based on this, the TITAN Penning trap is able to perform accurate 
mass measurements at a level of precision of below one ppb. 
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